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ABSTRACT 

The X-ray emission from most accreting white dwarfs (WDs) in symbiotic binary stars is quite soft. 
Several symbiotic WDs, however, produce strong X-ray emission at energies greater than ~ 20 keV. 
The Swift BAT instrument has detected hard X-ray emission from 4 such accreting WDs in symbiotic 
stars: RT Cru, T CrB, CD -57 3057, and CH Cyg. In one case (RT Cru), Swift detected X-rays out 
to greater than 50 keV at a > 5er confidence level. Combining data from the XRT and BAT detectors, 
we find that the 0.3-150 keV spectra of RT Cru, T CrB, and CD -57 3057 are well described by 
emission from a single-temperature, optically thin thermal plasma, plus an unresolved 6.4-6.9 keV 
Fe line complex. The X-ray spectrum of CH Cyg contains an additional bright soft component. For 
all 4 systems, the spectra suffer high levels of absorption from material that both fully and partially 
covers the source of hard X-rays. The XRT data did not show any of the rapid, periodic variations 
that one would expect if the X-ray emission were due to accretion onto a rotating, highly magnetized 
WD. The X-rays were thus more likely from the accretion-disk boundary layer around a massive, 
non-magnetic WD in each binary. The X-ray emission from RT Cru varied on timescales of a few 
days. This variability is consistent with being due to changes in the absorber that partially covers 
the source, suggesting localized absorption from a clumpy medium moving into the line of sight. The 
X-ray emission from CD —57 3057 and T CrB also varied during the 9 months of Swift observations, 
in a manner that was also consistent with variable absorption. 

Subject headings: stars: symbiotic — stars: individual (RT Cru, CH Cyg, CH Cyg, CD —57 3057) 



1. INTRODUCTION 

To an optical astronomer, a symbiotic star is a red 
giant with an additional hot blue spectral component 
and emission lines. In most cases, the blue component 
is produced by an accreting white dwarf (WD), making 
symbiotic stars cousins of cataclysmic variables (CVs; 
in which a WD accretes from a Roche-lobe filling, late- 
type dwarf). In the soft X-ray regime, many symbi- 
otic s tar s have been de t ected with ROSAT (jMiirset et alJ 
I1997D - IMiirset et all (|1997f) found that some symbi- 
otics could be classified as super-soft X-ray sources, in 
which the spectrum is approximately that of a black- 
body with kT ~ 50 eV and the emission emanates from 
the phot osphere of a hot WD with surface nuclear shell 
bur ning jK ahabka fe van den Heuvell I1997T ). For oth- 
ers, IMiirset et alJ (|1997D fitted the spectra with optically 
thin thermal plasma models and found temperatures of 
~ 10 7 K. They suggested that the emitting plasma was 
heated by the collision of the slow, dense wind from the 
M giant and the faster (^1000 km s _1 ) wind from the 
accreting WD. 

In this paper, we report on Swift (Gehrels et al. 
I2004D observations of 4 WD symbiotics that have been 
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detected by the Swif t Burst Alert Telescope (BAT, 
iBarthelmv eTall [2005) at energies above ~ 20 keV: 
RT Cru, T CrB, CD -57 3057 (sometimes referred to 
as SS73 17), and CH Cyg. Although we have pub- 
lished preliminary reports on the detection of these BAT 
sources and the subsequent identification with symbi- 
otic stars, here we present the first spectral and tem- 
poral analysis of the BAT data on these symbiotic stars. 
These observations reveal a new type of X-ray emission 
from WD symbiotics, since the extrapolation of super- 
soft and colliding-wind emission falls far short of the 
Swift BAT detections. After presenting the BAT survey 
data in jQ] and the p ointed Swift X-Ray Telescope (XRT, 
iBurrows et al. 2005) observations and data analysis in ^3] 
and fJH respectively, we summarize the combined XRT 
and BAT results in Sj5j In fj6l we consider the origin of 
the hard X-rays in these 4 symbiotic stars through com- 
parison with hard X-ray emission from CVs. 

2. SYMBIOTIC STAR IDENTIFICATIONS AND BAT 
RESULTS 

All BAT data in this p aper derive from the 9-m onth 
BAT survey described by IMarkwardt et all (|2005Q . In 
this survey, BAT data were collected with a nominal 
5-minute sampling frequency. The data retained their 
full spatial and energy resolution. The data were fil- 
tered to remove any contamination from high detector 
background rates due to South Atlantic Anomaly radi- 
ation, bad spacecraft pointing, and electronics noise in 
the BAT detectors. To flat-field the images, the con- 
tributions from known bright sources and smooth vari- 
ations due to sky background were subtracted. Sky im- 
ages were t hen produced using a cross-correlation tech- 
nique (e.g. iSkinner et al.|[l987l) . Finally, individual im- 
ages were combined into an all-sky mosaic, weighted by 
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the image noise level. These mosaics were then searched 
for excesses in each energy band, and a catalog of sources 
produced. To accurately identify BAT survey sources, a 
pointed observation with the Swift narrow field instru- 
ments was typically also obtained. Four WD symbiotics 
were detected in the BAT survey; we list their positions 
and count rates in Table??. 

RT Cru was detected in hard X-ra ys by the INTE- 
GRAL mission as I GR J12349-6434 (IBird et all [20071 : 
iMasetti et aT1l2005fl . iTueller et all (j2005d) confirmed the 
source identification with a Swift /XRT observation of the 
BAT survey source SWIFT J1234.7-6433. The pres- 
ence of various optical emission lines superimposed on 
an M giant continuum, as well as a UV excess, indicate s 
that RT Cru is a sym biotic star (|Cieslinski et al.|[l994T ). 
ICieslinski et all (j!994l ) further noted 0.03-0.05 mag opti- 
cal flickering with time scales of 10-30 minutes in the V 
band, a characteristic that RT Cru shares with a subset 
of symbiotic stars, including CH Cyg, T CrB, Mira AB 
(o Ceti), MWC 560 and RS Oph (ISokoloski et alJl200l . 

T CrB was detected in the BAT survey as 
SWIFT J1559. 5+2553. Subsequent observations with 
the Swift/XRT confirmed that this hard X-ray source 
was associated wi th the symbiotic star T CrB 
(|Tueiler et al.ll2005allbl ). T CrB is a recurrent nova (with 
eruptions in 1866 , 194 6) and has an orbital period of 
227.6 days ([Kraft 1958) . For a recen t study of its rapid 
variability, see IZamanov et aTl (|2004D . 

CD -57 3057 was detected in the Swift /BAT s ur- 
vey as SWIFT J1010.1-5747 (ITueller et alJ I2005allbh . 
It was also detected by the IBIS i nstrument on IN - 
TEGRAL (iReynivtsev et al.l l200l IBird et all 12001 . 
IMasetti et all (|2006t) suggested CD -57 3057 as the 
likely optical counterpar t of t he hard X-ray source. 
ISandu lcak & Stephenso n (1973) made the original de- 
ter mination that CD — 57 3057 is a symbiotic star, 
and|P ereira et al.l (|2003l ) performed additional studies. 
Smit h et all (|2008l ) reported on X-ray observations of 
this source with Suzaku. 

CH Cyg is the only object of the 4 that has been well 
studi e d in the trad i tional X-ray band (e.g. lEzuka et al.l 
T99l IMukai et all 120071 IGallowav fc Sokoloskil 12001 
Karovska et alJ 120071 ). With a distance of 245±40 pc 
(from the Hipparcos parallax; iPerrvman et al.lll997f ). it 
is one of the nearest symbiotic stars. Of the two known 
photometric periods, the 756-day period is thought t o 
be due to pulsation of the M giant (jSchmidt et al.1 12006) , 
and the 14.5-year period is probably the orbital period. 
CH Cyg was detected in the BAT survey. 

The BAT detected all 4 sources at energies above 
24 keV, with RT Cru detected above 50 keV with a con- 
fidence level of greater than 5a. Fitting the 4 channel 
BAT survey spectra to a simple bremsstrahlung model 
suggests that this hard emission is due to the presence of 
high-temperature plasma (kT ~ 15-30 keV for RT Cru, 
T CrB, and CD -57 3057, and - 5 keV for CH Cyg). 
TableQ] lists parameters of the fits to the BAT spectra 
for each source. Although the BAT spectra are not of 
high enough quality to make a strong statement about 
the hard X-ray emission mechanism, the lower \ 2 values 
for the bremsstrahlung model fits compared to a power- 
law model suggest that the bremsstrahlung model is pre- 
ferred. 

Each source varied during the 9-month BAT survey. 
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Fig. 1. — BAT light curves for 4 symbiotic stars during their Swift 
visibility period during 2005. Each source shows a clear variability 
in the hard X-ray band. Of particular note here is the decrease in 
the hard X-ray flux in CH Cyg, which was seen to be in a lo w flux 
state in a Suzaku observation reported by Mukai et al. (2007). 

Figure[T] shows the BAT survey light curves of these 4 
symbiotic stars during 2005, binned with approximately 
30-day time bins. Although there are no obvious pe- 
riodicities in the BAT light curves, fits to a constant 
model produced large values of reduced x 2 (x 2 , where 
v is the number of degrees of freedom; in this case xt 
ranged from 3.5 to 8.0 and v = 9). The hard X-ray 
flux from CH Cyg shows a clear decay to almost zero, 
which is consistent with th e low flux seen in 2006 with 
Suzaku (|Mukai et al.l[2007h . In 2000, INTEGRAL mea- 
sured a 10-50 keV flux of 3.4xl 0~ 10 ergs cm _2 s _1 for 
RT Cru (jChernvakova et ahll2005l ). which is substantially 
brighter than seen in the BAT measurements, indicating 
a dramatic X-ray variability over the ~5 years previous 
to the BAT measurements. 

3. POINTED Swift OBSERVATIONS 

As part of a program to use Swift narrow-field instru- 
ments to identify the counterparts of new hard X-ray 
sourc es from the BAT All Sky Survey (fMarkward t et all 
2005), Swift performed pointed observations of T CrB, 
CH Cyg, and CD -57 3057. RT Cru was initially ob- 
served to determine whether it was the X-ray counter- 
part to INTEGRAL source IGR J12349-6434. After the 
soft X-ray and optical counterparts to each of the 4 new 
hard X-ray source had been identified, Swift performed 
additional observations to obtain more detailed timing 
and spectral information. A log of these observations 
are presented in Tabled 

For all pointed observations, the XRT was in "Auto" 
state, in which the data collection mode is automati- 
cally chosen, based on the brightness of the source, to 
minimize the effects of pile-up. The two main modes of 
observing with the XRT are "Photon Counting" (PC ) 
and "Window Timing" (WT) modes (jHill et al.ll2005T) . 
PC mode provides full-field imaging and spectroscopy, 
with low temporal resolution (2.5 s). WT mode pro- 
vides higher temporal resolution (1.8 ms) at the expense 
of imaging. The RT Cru, T CrB, and CD -57 3057 
observations were taken exclusively in PC mode. For 
CH Cyg, the one pointed observation was performed in 
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TABLE 1 

Parameters of bremsstrahlung and power law model fits to the BAT 9-month survey 4-channel spectra 



Source 


kT 
(keV) 


Xt (dof) 


Photon Index 


Xt (dof) 


15-150 keV Flux 

(erg/s/cm 2 ) 


RT Cru 




0.27(2) 


9 o+U.3 
z -°-0.2 


2.56(2) 


5.6 x lO" 11 


T CrB 


17-2±|;1 


0.48(2) 


n 7+0.4 
z ''-0.3 


1.42(2) 


5.7 x lO" 11 


CH Cyg 


4.911:? 


0.80(2) 


4 fi +1S) 


0.88(2) 


2.6 x lO" 11 


CD -57 3057 


15.21^ 


0.93(2) 


9 q+0.6 


1.91(2) 


2.9 x lO" 11 



NOTE. — All errors are quoted at 90% confidence. 

a mixture of WT and PC modes. Since the XRT has only 
a thin optical blocking filter, however, the PC-mode data 
for CH Cyg were severely contaminated by optical load- 
ing on the XRT CCD from the bright M giant. When 
CH Cyg was observed in PC mode, the data from all but 
the wings of the XRT point spread function were there- 
fore effectively useless. The faster-exposure WT-mode 
data for CH Cyg were less affected by optical contami- 
nation and therefore suitable for spectral analysis. 

4. DATA ANALYSIS 

We used the Swift Release 2.5 software 8 to analyze the 
data, and we reduced the XRT data in the standard fash- 
ion with xrtpipeline. We determined the coordinates 
for XRT point sources with xrtcentroid, and performed 
all timing and spectral analyses with counts extracted 
from a circular region with a 20-pixel radius around the 
centroid of the XRT point source. A region with this ra- 
dius contains approximately 99% of all counts expected 
from an XRT point source. To correct for the presence 
of dead columns on the XRT CCD during timing anal- 
ysis of XRT data, we used the standard tool xr tlcorr. 
For s pectral fitting, we used XSPEC vl2.3.0 (jArnaudl 
1996) and the version-9 XRT calibration products. We 
performed simultaneous fits to the spectra from the cu- 
mulative XRT observations and those from the 9-month 
BAT survey observations for each source. Combining the 
XRT and BAT data allowed us to fit spectra spanning 
0.3 to 150 keV. 

5. RESULTS 
5.1. Timing analysis 

To search for periodic brightness modulations arising 
from the WD spin, we performed timing analysis on 
barycentrically corrected XRT event data for all 4 sym- 
biotic stars. Since WD spin periods typically lie in the 
range from mi nutes to a few hours, we calculated the Z\ 
statistic (e.g.. iBuccheri et aT1ll983|) for each source for a 
range of periods from 1 minute to 10 hours. We did not 
find any significant periodicities in the light curves of any 
of the symbiotic stars in our sample. 

In order to determine the likelihood of detection of 
a WD spin period in our data, we performed Monte- 
Carlo simulations to recreate our data with varying de- 
grees of modulation and periodicities in the 100 to 5000s 
rang e, following the most typ ical range of WD spin peri- 
ods ([Norton fc Watson 1989). We found that amplitude 
modulations of > 20%, where the amplitude modulation 
is defined as the sinusoidal amplitude divided by the peak 
flux, were detected at high significance, utilizing a Z\ 

8 http:/ /swift. gsfc.nasa.gov/docs/software/lheasoft/ 



search, in greater than 99% of our simulated datasets. 
iNorton fc Watsonl (|1989| ) report that X -ray spin modu- 
lations in magnetic CV systems are typically in the 40- 
90% range, which would be clearly detected in our data. 
Therefore if WD spin modulation is present in any of 
these systems, it is at a level of modulation far lower 
than seen in well-established magnetic WD systems. 

The sparsity of the pointed XRT data (see TableEJ did 
not permit a search for modulations with long periods 
(e.g., periods on the order of days or months). As the 
BAT light curves in Figure[T] show, however, the hard 
X-ray emission (15-25 keV) from all 4 symbiotic stars 
varied by a factor of 3 or more on a timescale of months. 
In addition, RT Cru showed considerable variability in 
the XRT band on a timescale of days. For example, 
between 2005 October 20 and 2005 October 22, the 0.3- 
10 keV flux tripled. On 2005 November 2, RT Cru was 
detected at 0.35 counts s _1 , which was a factor of ~6 
higher than the lowest flux level, on 2005 October 20. 
The 0.3-10 keV flux from CD -57 3057 also varied over 
the course of the Swift observations, but by at most a 
factor of 2. Also in contrast to RT Cru, the 0.3-10 keV 
flux from T CrB varied by only a relatively modest factor 
of 1.4 (maximum/minimum). We consider the associated 
changes in the spectral shape in i )5.3l 

5.2. Spectral Analysis 

As a first approximation for the 0.3-150 keV spec- 
tra of these 4 WD symbiotic stars, we fit an absorbed, 
single-temperature thermal bremsstrahlung model, with 
the addition of a Gaussian line component at around 
6.4-6.9 keV. Although plasma models that include spec- 
tral lines are more likely to correctly describe the data 
than simple thermal bremsstrahlung models, the qual- 
ity of our data did not allow us to detect complex line 
features, and plasma models did not provide improved 
fits over bremsstrahlung models. We thus parameter- 
ized all spectral fits using absorbed single-temperature 
bremsstrahlung models. Such a model provided a good 
fit for T CrB. For CD -57 3057 and RT Cru, however, 
absorbed single-temperature models had problems be- 
low 2 keV. The single-temperature fit to CH Cyg gave a 
strong low temperature component (kT ~ 0.2 keV) and 
significant residuals at higher energies, confirming the 
suggestion from the BAT data that additional broadband 
spectral components are needed to describe the spectrum 
above 2 keV. 

5.2.1. RT Cru 

An absorbed, single-temperature bremsstrahlung 
model provided a poor fit to the spectrum of RT Cru 
(xt = 4.0, v = 81). Because the initial fit generated 
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TABLE 2 

Log OF Swift POINTED OBSERVATIONS of symbiotic stars utilized in this paper. 



Target 


Observation ID 


Date of Obs. 


Exposure 
(ks) 


Mean XRT Count rate 
(count/s, 0.3-10 keV) 


HT (Jru 


00035099001 


2005-08-20 


4.5 


0.236 ± 0.009 




00035099002 


2005-10-20 


1.9 


0.058 ± 0.007 




00035099003 


2005-10-22 


18.9 


0.161 ±0.037 




00035178001 


2005-11-02 


2.2 


0.203 ± 0.012 




00030840001 


2006-12-18 


2.7 


0.17 ±0.01 




00030840002 


2006-12-21 


1.6 


0.154 ±0.007 


T CrB 


00035171001 


2005-06-17 


8.8 


0.063 ± 0.034 




00035171002 


2005-10-07 


5.2 


0.069 ± 0.005 




00035171003 


2005-10-19 


10.3 


0.095 ± 0.004 


CH Cyg 


00035239001 


2005-11-03 


9.2 


0.295 ± 0.009 


CD -57 3057 


00035324001 


2005-08-15 


3.2 


0.048 ± 0.006 




00035324002 


2005-08-17 


3.5 


0.062 ± 0.007 




00035324003 


2005-09-27 


4.6 


0.104 ±0.006 




00035324004 


2005-11-22 


3.8 


0.075 ± 0.006 



a soft excess in the residuals, we tried including a sec- 
ond absorbed bremsstrahlung component. The addition 
of the second broadband emission component resulted 
in a much-improved fit (x 2 — 1-2, v — 79). In this 
second fit, however, the two bremsstrahlung tempera- 
tures were consistent to within the errors (kT ~ 38 
keV), with only the absorption differing significantly for 
the two bremsstrahlung components: A^hi = (7.6^ 2 'i) x 

I0 22 cm" 2 and N m = 0.7t°:iJ X 10 22 cm- 2 . These re- 
sults strongly suggest the presence of an absorber that 
only partially covers the source. 

A single-temperature thermal bremsstrahlung model 
with one absorber that fully covers the source, and an- 
other that only partially covers the source (i.e., the 
XSPEC pcf abs model), provided a good fit to the data 
(xf, = 1.17, v = 80), with a fitted partial covering ab- 
sorption is Nh{PC)—^ 7 x I0 22 cm -2 and an intrin- 
sic absorption of N-r(FC) =~ 7 x I0 21 cm- 2 , consistent 
with the expect ed total Galactic absorbi ng column of 
9.6 x I0 21 cm" 2 (|Dickev fe LockmanlH990l ). Tabled lists 
the parameters of this fit. The temperature from the fit 
to the combined XRT and BAT spectrum is consistent, 
to within the errors, with the temperature from the fit 
to the BAT spectrum alone. Allowing the normaliza- 
tion for the bremsstrahlung fit to differ for the XRT and 
BAT data did not improve the fit. Moreover, the nor- 
malizations for the XRT and BAT data were consistent 
to within the errors. 

The XRT spectrum of RT Cru included a strong emis- 
sion feature associated with the 6.4-6.9 keV Fe line com- 
plex. To parameterize this feature, we fit the 5-10 keV 
portion of the XRT spectrum with a single-temperature 
bremsstrahlung continuum plus a variety of Gaussian 
lines. Although a single broad Gaussian provided an ad- 
equate fit to the line complex (% 2 = 1.140 for v = 40), 
both the physical parameters (a line energy of 6.63^q'qs 
keV and a Gaussian a of 269^^ eV) and the residuals 
suggested that the line complex actually required mul- 
tiple Gaussian components. Fitting a 2-Gaussian model 
gave a much improved fit (% 2 = 0.75 for v — 39), in 
which the line complex was well described by a narrow 
(cr < 90 eV) Fe Ka line at 6.461^5 keV and a broader 
{a = 178tg|° eV) line, probably due to H-like Fe xxvi 

emission, at 6.911q \g keV. Given the detection of Fe XXV 
line emission at 6.7 keV in the Chandra grating spectra 
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Fig. 2. — Two component Gaussian model fit to the RT Cru Fe 
line complex. Dashed lines show the 3 components of the fit, a 
bremsstrahlung continuum model and two Caussian lines repre- 
senting a narrow 6.4 keV Fe Ka line and a broad 6.9 keV line of 
H-like Fe. 

of RT Cru (jLuna fe Sokoloskil[2007l ). in addition to the 
Fe Ka and Fe xxvi lines suggested by the XRT fits, it 
is likely that the width of the 6.9 keV line is due to 
blending of the 6.7 and 6.9 keV lines. A fit to the XRT 
spectru m with 3 line component s fixed at the energies 
given bv lLuna fe Sokolo ski ( 2007) provided a good fit to 
the data, but no statistical improvement over the two- 
line model. Figure[2] shows the two-line fit to the Fe-linc 
complex. We provide additional discussion of the spec- 
tral continuum, which is variable, in M5 - 31 

5.2.2. T CrB 

For T CrB, the single-temperature bremsstrahlung 
model with a Gaussian line provided an acceptable fit to 
the full 0.3-150 keV spectrum (xl = 1.1 with v = 30). 
T CrB has a high level of absorption, with N-r(FC) — 
(21.5^2 3) x 10 22 cm- 2 . The total Galactic absorbing 
column in the direction of T CrB is 4.7 x 10 20 cm -2 
(jDickev fe Lockmanl[l990l ]. The fitted value is more than 
400 times this maximum Galactic value, and therefore 
the majority of Nn(FC) must be intrinsic to T CrB. 
The temperature from our fit to the full spectrum was 
kT = 28.0^5 '\ keV — somewhat higher than the tem- 
perature from the BAT spectrum alone (Tablc[T]). The 
Gaussian line in this model is required with a high signif- 
icance; the fitted line energy is 6.6 ±0.1 keV, with a line 
width of a = 400±^ eV. The large width of this line 
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TABLE 3 

Partial covering model fits to combined BAT and XRT spectra. 

Source N H (FC) a N H {PC)» PCF C kT Line Peak Line Width xt (dof) 

xlO 22 cm 2 xl0 22 cm 2 (kcV) ( kcV ) . ( cV ) 

RT Cru 0.73^;^ 6.90±^ 0.83^;^ 37.lj£| 6.41± ;£> 6.06^^;^ 1.17 (80) 

T CrB 0.047(fixed) 27.8±^ 0.96±°;°f 17-2j;f;° 6.64±°;^ 421± 2 °* 0.84(28) 

CD -57 3057 2.03±g;|g 31.8±j;|';| 0.87±°;°^ 17.3±|;g 6.61±g;^| 180^ 1.34(43) 



NOTE. — All errors are quoted at 90% confidence. 

a ./Vf/(FC) is the column density of absorbing material that fully covers the source. 
JVjj(PC) is the column density of absorbing material that partially covers the source. 
C PCF is the partial-covering fraction. 
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Fig. 3. — Absorbed bremsstrahlung plus Gaussian line model fit 
to the combined XRT (black) and BAT (red) data for T CrB. Note 
that the absorption is a combination of the line of sight Galactic 
absorption and a partial covering fraction model to account for 
intrinsic absorption. 

suggests the presence of a number of blend ed Fe lines, 
like those seen in CH Cyg (jMukai et al.ll2007t) . The XRT 
spectrum for this source, however, did not have a high 
enough signal-to-noise ratio to identify multiple lines. 

As the XRT and BAT observations were not contem- 
poraneous, and T CrB varies in both the XRT and 
BAT energy bands, we re-fit the 0.3-150 keV spectrum 
with the normalizations for the XRT and BAT data free 
to differ. With all other parameters forced to be the 
same for the two data sets, we obtained an improved fit 
(xt — 0.94 for v = 29) with a bremsstrahlung temper- 
ature of kT = I7.6I5 3 keV, which was consistent with 
the temperature from the fit to the BAT spectrum alone. 
The other parameters remained unchanged. The normal- 
ization for the BAT data was approximately 75% higher 
than that for the XRT data, suggestive of variations in 
brightness between the time of the BAT and XRT obser- 
vations. Allowing kT — as opposed to the normalization 
- to have different values for the XRT and BAT data 
sets did not provide a significant improvement over sim- 
ply allowing the BAT normalization to vary. There was 
therefore no evidence for spectral variations within the 
BAT energy band. 

We also fit the spectrum with a model that included 
partial-covering absorption, with the column density of 
the absorber that fully covered the sourced fixed to the 
Galactic value of 4.7 x 10 20 cm~ 2 (we note however, that 
the actual line of sight absorption is most likely lower 
than this) . The resulting fit was an improvement over the 
absorbed bremsstrahlung model (xt = 0.84 for v = 28) 
with a fitted partial covering absorption is A#(PC)=^ 
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Fig. 4. — Absorbed bremsstrahlung plus Gaussian line model fit 
to the combined XRT (black) and BAT (red) data for CD -57 3057. 
Note that the absorption is a combination of line of sight absorp- 
tion and a partial covering fraction model to account for intrinsic 
absorption. 

3 x 10 23 cm~ 2 . An F-test comparing this model with 
the full-covering absorber model gave a 4% probability 
of this improvement being seen in a random data set. 
Thus, although the fitted covering fraction was close to 
100%, the partial-covering absorber model was still the 
better one for T CrB. We present the results of this fit 
in Figure[3] and Tablegj 

5.2.3. CD - 57 3057 

Fitting the combined XRT and BAT spectra of 
CD —57 3057 with a model that consisted of a single- 
temperature thermal bremsstrahlung plus a Gaussian 
line, with partial-covering absorption, gave an accept- 
able fit (xl = 1-34 for v — 43). Figured] shows the data, 
model spectrum, and residuals, and Table[3] lists the pa- 
rameters of the fit. The column density of absorbing ma- 
terial that fully covered the source was high (iV/j-(FC)=~ 
2 x 10 22 cm -2 ), but consistent with the total Galactic ab- 
sorption of 1.7 x 10 22 cm -2 ()Dickev fe Lockmanlll990l) . 
The column density of absorbing material that only par- 
tially covered the source, iVg(PC) =~ 3 x 10 23 cm" 2 
was similar to the value for T CrB. ISmith et all (|2008l ) 
reported detections of Fe Ka, Fe xxv and Fe xxvi lines 
in the Suzaku spectrum of CD —57 3057; the emission- 
line energy and width that we measured in the Swift, 
XRT spectrum were therefore most likely the result of 
the blending of these lines. 

5.2.4. CH Cyg 

As there was only one pointed observation of CH Cyg 
during 2005, the signal-to- noise ratio for the resulting 
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Fig. 5. — Hardness ratio as a function of 0.3-10 keV brightness 
for RT Cru, T CrB, and CD -57 3057. The hardness ratio is the 
ratio of the 4-10 keV flux and the 0.3 - 4 keV flux. Note the 
strong anti-correlation between X-ray brightness and the hardness 
of the source. 

Swift XRT spectrum was low. Moreover, due to the 
optical brightness of this object, a large portion of the 
observation was performed in WT mode, further reduc- 
ing the effective exposure time. Nonetheless, the XRT 
spectru m appeared consiste nt with that found by Suzaku 
in 2006 (Mukai et al. 2007), in which the spectrum was 
dominated by a soft (E < 2 keV), relatively unabsorbed 
component, but also contained an additional hard (E > 
2 keV), highly absorbed component (Njj > 10 23 cm -2 ). 

5.3. Time dependent spectral variability 

To investigate the spectral variability of the sources 
with multiple observations (i.e., RT Cru, T CrB, and 
CD —57 3057), we used the XRT data to calculate hard- 
ness ratios as a function of time. For these hardness ra- 
tios, we selected the bands 0.3-4 keV and 4-10 keV. Divi- 
sion of the XRT spectra at 4 keV gave an approximately 
equal number of photons in each band for RT Cru. For 
consistency, we used the same energy bands for T CrB 
and CD -57 3057. 

RT Cru varied significantly on a timescale of days, in 
one case showing a variation of a factor of 3 in observa- 
tions over a 5-day period. Moreover, the hardness ratio 
and X-ray brightness appear to have an inverse relation- 
ship. Figure[5] shows the hardness ratio as a function of 
brightness for RT Cru, T CrB, and CD -57 3057. For 
RT Cru and CD —57 3057, there is a clear inverse re- 
lationship between hardness ratio and brightness. For 
RT Cru, the linear correlation coefficient is —0.85 (for 
15 data bins). The relationship between hardness and 
XRT count rate for RT Cru is, however, actually best 
described by a power-law with an index of —1.2, so the 
correlation is most likely better than the linear coeffi- 
cient suggests. Although T CrB also shows a suggestion 
of this same inverse relation, the relatively low dynamic 
range of X-ray brightness measured for this object leads 
to a correlation that is not statistically significant. 

The nature of the variations from RT Cru and 
CD —57 3057 suggests that they could have been due 
to changes in the spectral shape of the soft portion of 
the X-ray spectrum, rather than changes in the over- 
all normalization. To determine the nature of the ob- 
served hardness-brightness correlation, we divided the 
XRT data for RT Cru into three brightness ranges and 



Fig. 6.— Model fits to combined XRT and BAT RT Cru spec- 
trum at 2 flux levels. Black points are XRT data of the source 
in its high state, red points the XRT spectrum at intermediate 
brightness, green points are BAT data. A simple brcmsstrahlung 
plus Gaussian line model with both intrinsic and partial covering 
absorption fits these data. For the intermediate and high states, 
only the partial covering fraction is allowed to vary in the spectral 
model. The normalization of the BAT data is independent from 
the XRT to compensate for instrumental differences and systematic 
differences due to time coverage. 

extracted spectra for each brightness. Since the RT Cru 
light curve at 8-10 keV was essentially constant, the most 
obvious way to generate changes in the spectra shape was 
variable absorption. We thus fit the spectra for the dif- 
ferent brightness states simultaneously using the partial- 
covering model, with all the parameters fixed to the 
values from Table[3] except for the partial-covering frac- 
tion, which we left as a free parameter. This procedure 
produced good fits for the spectra from each brightness 
state, with the fitted partial covering fractions values of 
0.64t°;?4, 0.95tg;$i, and a 90% confidence lower limit of 
0.99 for the high, intermediate, and low brightness states 
respectively. FigureJS] shows the spectra corresponding to 
the intermediate and high states. 

The simplest way to explain the variations in RT Cru 
is with a variable partial-covering fraction, perhaps due 
to clumpy absorbing material moving in and out of our 
line of sight. We suggest that such a process causes the 
brightness variations in RT Cru, CD —57 3057, and pos- 
sibly also in T CrB, although the significance of the spec- 
tral variability in the latter two sources makes it less 
clear, especially for T CrB. However, their spectral simi- 
larities, and the nature of their variability, strongly sug- 
gests that RT Cru, CD -57 3057, and T CrB belong to 
a subclass of highly absorbed, hard X-ray emitting sym- 
biotic stars that do not show any evidence for a CH Cyg- 
like soft component. 

6. DISCUSSION 

Previous X-ray observations of WD symbiotic stars re- 
vealed primarily soft X-ray emission. Utilizing the hard 
X-ray imaging capabilities of INTEGRAL and the Swift 
BAT, a new class of hard X-ray emitting WD symbiotics 
has been identified, of which RT Cru, CD -57 3057, 
T CrB and CH Cyg are members. 

The hard X-ray component appears to be predomi- 
nantly thermal in origin, judging by the bremsstrahlung 
continuum shape (BAT; this work) and the presence 
of Fe K lines. Moreover, deeper observations re- 
ported elsewhere (|Ezuka et al.l 119981 : iMukai et al.l 120071 ; 
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iLuna fc So koloski 2007) establish the presence of the Ka 
lines from H- and He-like Fe, and our XRT data are 
consistent with the presence of these lines. The hard 
X-ray component of all 4 objects is strongly absorbed, 
and an absorber that only partially covers the source is 
required in all but T CrB. We do not detect any peri- 
odic variations. We did, however, detected strong aperi- 
odic variability in RT Cru, CD —57 3057, and perhaps 
T CrB, which spectral investigation showed to be caused 
predominantly by variability in the partial covering ab- 
sorber. 

We can estimate the 15-150 keV luminosities of these 
symbiotic stars from the BAT fit results (TableQJ if the 
distances are known. For CH Cyg, we use the Hipparcos- 
measured distance of 245 pc to infer a luminosity of 
~ 2 x 10 32 ergss -1 . Note, however, that during the 
period of the BAT survey observation the source bright- 
ness decreased to a hard X-ray low state seen by Suzaku 
(|Mukai et al.ll2007h . The distance to T CrB is estimated 
to be about 1 kpc based on its K magnitude, where in- 
terstellar extinction and the contribution fro m the accre- 
tion d isk is less important than in the optical (jHric et al.l 
1998). Based on this distance, the 15-150 keV luminosity 
of T CrB is - 7 x 10 33 ergss" 1 . 

For the other two systems, the distances are highly 
uncertain. The inferred luminosities are ~ 3.5 x 
10 33 [d/(l£;pc)] 2 ergss" 1 and - 7 x I0 33 [d/(lkpc)] 2 
ergss -1 , respectively, for CD —57 3057 and RT Cru. 
Given that CD -57 3057 and T CrB have similar 
V mag nitudes, 1 kpc is a reasonable guess for the 
former (|Smith et al.l 120081 . argue that the distance to 
CD —57 3057 is between 250pc and 1 kpc), although 
the accretion disk contributions and the interstellar ex- 
tinctions may we l l be different for these two stars. 
ILuna fc Sokoloskil (|2007D suggest a distance of 1.5-2.0 
kpc for RT Cru, which would mean that RT Cru 
has an intrinsic luminosity 4 to 8 times brighter than 
CD -57 3057, or greater if CD -57 3057 is closer than 
1 kpc. 

Considering these luminosity estimates, and the ther- 
mal nature of the X-ray spectra, we believe that the hard 
X-rays in these symbiotic stars are due to accretion onto 
WDs. The luminosity of such a component is propor- 
tional to the gravitational potential of the white dwarf 
(GM/R, where G is the Newtonian constant of gravity, 
M is the white dwarf mass, and R its radius) times the 
mass accretion rate. The highest temperature in the ac- 
creting plasma is also primarily controlled by GM/R, al- 
though it also depends on the de tails of the accretion ge - 
ometry. Based on equation (5) of lLivio fc Warner! (1984), 
these symbiotic stars could have large (~ 10 12 — 10 13 
cm) disks. Although such large disks differ markedly 
from CV disks on AU scales, on small scales where ac- 
cretion energy is expected to be radiated as X-rays, the 
astrophysics is similar to that of CVs, and the physics 
of accretion is governed by the same fundamental laws. 
We therefore proceed by applying the knowledge gained 
from CVs to symbiotic stars. 

On scales of 10 9 cm, the magnetic field of the white 
dwarf can play a crucial role in shaping the X-ray emis- 
sion (jMukai 20051). Accretion in magnetic CVs proceeds 
vertically, forming a strong shock. Accretion in non- 
magnetic CVs proceeds via a Keplerian disk. The in- 



terface between the disk and the WD ("the boundary 
layer" ) is thought to be the origin of the X-rays in non- 
magnetic CVs. The X-ray luminosities of non-magnetic 
CVs are generally in the range 10 30 -10 32 ergss - typ- 
ically with kT~5- 10 keV bremsstrahhmg-like spectra 
(Bask ill et a,l.ll2005h . Magnetic CVs are often harder and 
more X-ray luminous: iSazonov et alj (|2006l ) derived lu- 
minosities in the range 3 x 10 31 -3 x 10 33 ergs s -1 for mag- 
netic CVs detected by the RXTE slew survey. The 4 WD 
symbiotics that we detected with BAT appear to have 
luminosities that are comparable to, or slightly greater 
than, those of magnetic or non-magnetic CVs. 

The magnetic nature of a CV can be established by 
the detection of the spin period in the X-rays and/or the 
optical. That we did not find any coherent periodic vari- 
ations in the Swift XRT data, and that none have been 
reported in other data sets, is a strong indication against 
the magnetic nature of the accreting WDs in these sys- 
tems. That non-magnetic CVs are generally fainter and 
softer than magnetic CVs raises the question of why these 
symbiotic stars are luminous enough and hard enough to 
be detectable with INTEGRAL and Swift BAT. We note 
that only one non-magnetic CV has been detected by IN- 
TEGRAL: SS Cyg, which is the brightest non -magnetic 
CV in the 2-10 keV band (iBarlow et al.ll20M ). 

We argue that the high temperature and high luminos- 
ity of the 4 symbiotic stars are most readily explained 
by high values of GM/R, i.e., high mass white dwarfs. 
In the following, we estimate the masses of the WDs in 
T CrB, CD -57 3057, and RT Cru by assuming that 
the plasma temperature in a boundary layer around a 
WD of a certain mass is related to the plasma temper- 
ature in a magnetic accretion column onto a WD of the 
same mass by one of two simple scale factors. We obtain 
empirical estimates of these scale factors by comparing 
the measured bremsstrahlung temperatures in SS Cyg, 
a non-magnetic CV with an independently determined 
mass estimate, to the theoretical shock temperature for 
a magnetic CV of that mass. The vertical accretion in 
magnetic CVs allows plasma to reach a shock temper- 
ature of kT Stm = 3GM firriH /&R, where kT s ^ m is the 
shock temperature, m# the mass of a hydrogen atom, 
and n the mean molecular weight (we adopt 0.615 ap- 
propriate for Solar abundance plasma). When one fits 
a single-temperature bremsstrahlung model, the result- 
ing temperature is lower than kT s ^ m because the emis- 
sion is actually from a multi-temperature plasma cooling 
from kT S:Tn to ~0 keV. Nevertheless, single-temperature 
bremsstrahlung fits returns values in the kT ~ 1 0-50 
keV range for magnetic CVs (see, e.g., IIshidallT99ll for a 
compilation of Ginga results) . 

The shock temperatures in non-magnetic CVs are 
lower than in magnetic CVs for two reasons. First, 
half of the gravitational potential energy has been ra- 
diated away by the accretion disk. Second, the shocks 
in the boundary layer are oblique, and not as strong 
as those in magnetic CVs. Unfortunately, the details 
of the accretion disk boundary layers are complex and 
theories are not yet sufficiently advanced to provide 
an unambiguous guidance. Empirically, SS Cyg has a 
bremsstrahlung te mperature of 17 keV i n quiescence and 
7 keV in outburst (Whcatley et al. 2003) associated with 
the transition of the boundary layer from optically thin 
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to pr edominantly optically thick ((Patterson fc Raymond! 
11985ft. Its WD mass is thou ght to be M = 1.1 ± 0.1 M 
( Mauche fc Robinsonl |20Q~TI and references therein), and 
therefore its radius is 4900+/-700 km using the standard 
mass-radius relati onship (we have made use of the ana- 
lytical formula of IPringle fc Webbinkl [T975h . and there- 
fore kT s rn — 71 keV. Let us proceed, for now, by assum- 
ing that there is one scaling factor for the temperature 
obtained through a bremsstrahlung fit for low accretion 
rate (quiescence) boundary layer, and another for high 
accretion rate (outburst) boundary layer. These factors 
therefore are 17/71=0.24 for quiescence and 7/71=0.10 
for outburst. 

T CrB and CD —57 3057 have a bremsstrahlung tem- 
perature near 17 keV. Using the scaling factors above, 
the inferred fcT s , m is 17/0.24=71 keV or 17/0.10=170 
keV, depending on whether the boundary layer in these 
symbiotics resemble that of SS Cyg in quiescence or in 
outburst. These values imply the WD mass of 1.1 M© or 
1.35 M© for these two cases, respectively. Similarly, the 
30 keV bremsstrahlung temperature of RT Cru requires 
a 1.3 M or a 1.42 M© WD. The WD in CH Cyg, with 
its much lower temperature, is presumably less massive. 

Our estimate of the mass of the WD in T CrB in the 
high accretion rate case (1.35 M©) is consistent with 
those based on the recurrence time sca le of this recurrent 
nova (see, e.g.. H achisu fc Katolll999fl. Since the UV lu - 
minosity of T CrB is estimated by ISelvelli et al.1 (|1992l ) 
to be 40 L© or 1.6 x 10 35 ergss -1 , with an inferred mass 
accretion rate of ~ 2.5 x 10 -8 M© year -1 , it indeed must 
be in the high state. As in SS Cyg in out burst and in 
the old nova V603 Aql (|Mukai fc Oriol[2005l) . the UV lu- 
minosity of T CrB far exceeds its hard X-ray luminosity, 
probably because the bulk of the boundary layer is opti- 
cally thic k and predominantly emits in the UV and soft 
X-rays ()Patterson fc Raymond! [1985). Given the intrin- 
sic as well as interstellar absorption, the lack of detection 
of the soft X-ray component in T CrB is not surprising. 
Thus, scaling using SS Cyg in outburst leads to a solution 
that is consistent with previous studies of this recurrent 
nova system. 

All 4 systems show strong intrinsic absor ption. In ad- 
dition to the spectral results present here, lEzuka et al.l 
(1998) report that a partial covering absorber model is 
required to describe the strong intrinsic absorption seen 
by ASCA in CH Cyg. The variability of the intrinsic ab- 
sorption seen in RT Cru and CD —57 3057 over relatively 
short time scales (days) would suggest that the absorb- 
ing material is associated with accretion onto the WD, 
and cannot be explained by the mere presence of the M 
giant wind. Assuming Keplerian motion, a time scale of 
days corresponds to a distance from the WD on the or- 
der of 10 11 cm. These quantities are much smaller than 
the orbital time scale of years, and the expected binary 
separation, respectively. On the other hand, 10 11 cm is a 
reasonable distance for a feature associated with the ac- 
cretion disk. If the accretion disk is geometrically thin, 
we expect significant intrinsic absorption only in high 
inclination systems. For a thin disk, there is nothing ob- 



vious that would prevent us from observing hard X-rays 
from low inclination symbiotic stars, with little intrinsic 
absorption (absorption in the wind of the M giant is much 
less t han is observed in our 4 targets; Ivan den Berg et all 
2006). If the accretion disks in these high mass, high ac- 
cretion rate, symbiotic stars are geometrically thick, it 
may be possible for the majority of the sight lines to 
pass through local absorbers. 

We note that in the models of s tellar wind dis r uptio n 
by a compact object reported bv iBlondin et alJ <| 19901 ) . 
localized instabilities in the accretion flow can lead to 
absorption variability with timescales on the order of 
hours, consistent with timescales seen in our observa- 
tions. However these models also predict a strong vari- 
ability in absorption at the orbital period not seen in our 
data. As the orbital parameters of these systems are not 
well understood, comparison with such models is consid- 
ered speculative. 

7. CONCLUSIONS 

Swift BAT survey and pointed XRT observations have 
established CH Cyg, T CrB, RT Cru, and CD -57 3057 
as bright hard X-ray symbiotic stars. We have put forth 
an interpretation in which the hard X-rays are produced 
by accretion onto a massive, non-magnetic WD. Using a 
scaling factor from SS Cyg in outburst, we infer a WD 
mass of ~ 1.35 M© for T CrB. The same conclusion has 
already been reached by other, independent means. A 
1.35 M© white dwarf (i?=2700 km) has a gravitational 
potential well that is >7 times deeper than that of a 
0.6 M© (i?=8800 km) white dwarf. Other things being 
equal, the X-ray luminosity and temperature of such a 
massive white dwarf are therefore > 7 higher than those 
of a typical (0.6 M©) accreting white dwarf. The high 
mass of the WD in T CrB is thus the main reason that 
it is detected by Swift BAT. If we then take the hard 
X-ray emission from T CrB, which is known to contain 
a high-mass WD, as validation of our method, we can 
also infer near Chandrasekhar masses for the WDs in 
CD —57 3057 and in RT Cru. If these mass and accretion 
rate estimates are confirmed, these symbiotic stars are 
also likely to be recurrent novae. Moreover, given their 
near Chandrasekhar masses, they may be valid candidate 
progenitors of Type la supernovae. 

Higher sensitivity surveys are necessary to discover the 
full population of hard X-ray emitting WD symbiotics. If 
high absorption is universal, but the high temperature is 
not (as suggested by CH Cyg), then a wide-area survey 
in the 2-10 keV band may be necessary. At the same 
time, T CrB-like sources may stand out easily in deeper 
surveys of the >10 keV sky, including the future versions 
of the Swift BAT all-sky survey. 
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